Signature of effective mass in crackling noise asymmetry 
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Crackling noise is a common feature in many dynamic systems 
0, Q , the most familiar instance of which is the sound made by a sheet of paper 
when crumpled into a ball. Although seemingly random, this noise contains 
fundamental information about the properties of the system in which it occurs. 
One potential source of such information lies in the asym metric shape of noise 
pulses emitted by a diverse range of noisy systems llOUlll . 12] . but the cause 



of this asymmetry has lacked explanation [1]. Here we show that the leftward 
asymmetry observed in the Barkhausen effect Q - the noise generated by the 
jerky motion of domain walls as they interact with impurities in a soft magnet 
- is a direct consequence of a magnetic domain wall's negative effective mass. 
As well as providing a means of determining domain wall effective mass from 
a magnet's Barkhausen noise our work suggests an inertial explanation for the 
origin of avalanche asymmetries in crackling noise phenomena more generally. 

Crackling noise is the response of many physical systems to a slow external driving: 
outbursts of activity (avalanches or pulses) spanning a broad range of sizes, separated by 
quiescent intervals In condensed matter, notable examples are the magnetization noise 
emitted along the hysteresis loop in ferromagnets (i.e., the Barkhausen effect ]2j), the noise 
from magnetic vortices in type II superconductors |3j, ferroelectric materials 4] and driven 
ionic crystals ^]. In the context of mechanics, examples are the acoustic emission signal in 
fracture 



andplasticity []| and, on a larger scale, seismic activity in correspondence to an 
eartn^flQ. Quantitative understanding of cracklin g noise is of fundamental impor- 



tance in different applications, from non-destructive material testing to hazard prediction. 
This goal can be achieved only through the identification of general universal properties 
common to these systems, irrespective of their differences in the internal dynamics and mi- 
crostructural details. In this context, the average shape of the individual pulses composing 
the signal has been recently proposed as the best tool to characterize these universal fea- 

n 

tures of crackling noise jjj] . In analogy with critical phenomena, it is expected that pulses of 
different durations can be rescaled on a universal function, whose shape would only depend 
on general features of the physical process underlying the noise. This scenario is supported 
by the analysis of a variety of models, where pulse shapes are described by universal sym- 
metric scaling functions 3, 13. 11^. In most experimental data, however, the pulse shape 
is markedly asymmetric with respect to its midpoint, i.e. avalanches start fast but return 
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to zero more slowly llCt IllL Il2| . These results are puzzling because the models 

accurately reproduce several other universal quantities, such as avalanche distributions and 
power spectra [TlL [l5| . 

One of the most studied examples of crackling noise is the Barkhausen effect recorded in 
soft magnetic materials. The Barkhausen noise is due to the motion of domain walls, the 
interfaces separating regions of opposite magnetization, in response to a magnetic field (see 
the movie as Supplementary Information). Domain walls are characterized by an effective 



mass 



which is related to the increase of the wall energy with velocity, as experi- 



mentally revealed, for instance, in the dynamic susceptibility of insulating ferrite materials 
Q. far o m a gn e ts , _ inertia! effect s are usually neglected, being much 

smaller than eddy current dissipation [19(. This approximation is usually assumed in the 
description of the Barkhausen effect Q|. 

We show in this Letter that the asymmetry of pulse shapes is a direct signature of the 
effective mass associated to the objects moving under the action of the external field. In 
conducting ferromagnets the mass is negative and this results in a leftward asymmetry of 
the Barkhausen noise pulses. Our findings clarify the general meaning of pulse shapes in 
crackling noise phenomena: the form of avalanche displays an asymmetry that depends on 
the duration, and encodes important information on the characteristic time of the under- 
lying dynamics. Only on very long timescales pulse shape symmetry and universality are 
recovered. 

We measure Barkhausen noise pulse shapes in two ferromagnetic alloys. The experimental 
setup is depicted in Fig. ^ (see Methods section). In response to an increasing external field, 
domain walls move, inducing in the pickup coil a voltage signal v(t) proportional to the 
magnetization rate. The shape of the pulse is defined as the voltage (v(t,T)) at time t 
averaged over all avalanches of duration T. When appropriately normalized and plotted 
as a function of t/T, the experimentally measured shapes for different durations rescale 
fairly well (though not perfectly, see below) on a universal functionfsee Fig. El left), that 
is clearly asymmetric, in agreement with earlier measurements [1(1 llll. Il2|. Such a leftward 
asymmetry implies that pulses start rapidly and decay slowly. This is precisely the opposite 
of the effect of standard inertia, the resistance of a body to changes on its motion, which 
would imply a slow increase of the velocity when the domain wall is at rest. The asymmetry 
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can be quantified by computing the average skewness 

±S*dt{v(t,T))(t-i) 2 ] 

where t = 1/T J T dt(v(t,T))t. As shown in Fig. El in both samples the skewness is always 
positive, indicating a leftward asymmetry, and it displays a peak for Tp Xp pa 200/^s. 

To account for these experimental results, we make use of a simple and successful model 
for the Barkhausen effect in soft metallic ferromagnets, based on the dynamics of a planar 
domain wall in an effective pinning potential |2QJ]. The equation of motion for the wall 
position x is given by 

px = 2I s (c H t-kx + W(x)), (2) 

where (3 is a damping constant, I s is the saturation magnetization, cut is the external field 
increasing at rate cjj, —kx is the demagnetizing field and W(x) is a random field with 
Gaussian distribution and Brownian correlations j^. These correlations are believed to 
represent an effective description of a more general model with flexible domain walls (3, |^ . 
Eq. (J2J) can be solved exactly and provides an excellent description of the statistical properties 
of the Barkhausen noise, considering that the domain wall velocity x is proportional to the 
recorded voltage signal v(t). The solution of the model, however, yields a symmetric pulse 



shape 



13, 



141 ] . at odds with experimental evidence. 



The assumption used to derive Eq. (J2J) is that at each time the work done by the ef- 
fective field (i.e. the applied field corrected by the demagnetizing and pinning fields) is 
compensated by the e nerg y dissipated by eddy currents, which is estimated in the qua- 
sistatic approximation [20]. A more detailed analysis of eddy current dissipation, including 
dynamic effects (see the Methods section and Ref. t 23]), leads to the identification of a 
frequency-dependent effective mass, which turns out to be negative in the entire spectrum 
and equal to M* pa —[3t/(2j) at low frequencies, where r is the longest relaxation time 
and 7 ~ 1.05. For the materials considered in our experiments, the effective mass can be 
estimated to be M* pa —7 ■ 10~ 5 kg/m 2 , much larger than the positive Doring domain wall 
mass (M ~ 10" 9 'kg /m 2 This negative inertial effect can also be formulated by adding 

a non-local damping term in the equation of motion for the wall (see Methods section) 

Ti + — f e- (t - f ' )/T x{t')dt' = 2I s {c H t -kx + W(x)), (3) 
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where T and To are coefficients of the same order of magnitude with r + To = (3. In order to 
understand the role of the effective mass on the pulse shape, we have numerically integrated 
Eq. (jSJ) for different values of r. While the distributions of avalanches duration and size are 
unaffected by the addition of the inertial term, the pulse shapes become asymmetric and 
bear a remarkable similarity with the experimental ones (see Fig.|21 right). Also in this case 
the skewness is always positive (Fig. EJ), indicating a leftward asymmetry, and it displays a 
peak in correspondence with a characteristic time T p pa lOr. 

Using the value of the largest relaxation timescale for the samples considered experimen- 
tally (i.e. t = fiab 2 /7i 2 pa 5/xs, where b is the sample thickness, /i is the permeability and a 
the conductivity), we obtain a theoretical estimate of the peak position T p pa 50[is, which 
is reasonably close to the value measured experimentally, considering the approximations 
involved in the model. In particular, we have treated a single domain wall moving at the 
center of the sample, while in the experiments several domain walls are present. This in- 
duces a non trivial interference among the eddy current patterns generated by each wall. 
To corroborate our interpretation of the asymmetry in terms of a negative effective mass, 
we have performed simulations of Eq. (JHJ) replacing the non-local term with a positive mass 



term Mx 



2J| . The resulting skewness, reported in Fig. |1J is negative and displays a peak 



at T p pa 10M/r. 

In conclusion, we have shown that the pulse shape asymmetry commonly observed in 
Barkhausen noise measurements in metals is a signature of a negative effective mass of 
domain walls. Remarkably, the skewness exhibits a peak that can be used to track the 
characteristic relaxation timescale, corresponding to the ratio between mass and damping 
constant. Similar asymmetric pulse shapes are also observed in seismic data, recorded in 
correspondence to an earthquake j^, |^|. Seismic movements are due to the motion of fault 
planes in response to the stress accumulated in the crust. A non-local dynamical effect, 
analogous to the one discussed here, could be due to the presence of stress overshoots 25]. 
It would be interesting to test this mechanism using existing fault dynamics models. 



5 



Methods 



Materials and experimental methods 

In our experiments a long solenoid provides an homogeneous driving field H ramped at 
constant rate c# = 4A/(m s) , while a secondary pickup coil around the sample cross section 
measures the induced flux. The pickup coil is made of 50 isolated copper turns, wound within 
1 mm. Such a small width is required to avoid spurious effects due to demagnetizing fields. 
The measurements are performed only in the central part of the hysteresis loop around the 
coercive field, where domain wall motion is the relevant magnetization mechanism (1). In 
this case the recorded voltage signal v(t) is proportional to the domain wall velocity x. To 
reduce excess external noise during the measurement, we use a low pass pre- amplifier filter 
with a cutoff frequency of 100 kHz. We employ two ribbons of Fe64Co2iB 15 , having width 
a = lcm, thickness b = 24/im, and length c = 20cm: the first sample is amorphous, and 
subjected to a small tensile stress (2MPa). The second has been partially crystallized (5%) 
by thermal annealing. Tensile stress and thermal annealing are used to relax part of the 
internal stresses, leading to a high signal to noise ratio. 



Derivation of domain wall motion from eddy current dissipation 

Considering the configuration described in Fig. Q the relevant component of field H e = 
H e (x,y,t)z, associated with the eddy currents produced by the moving domain wall, obeys 
the equation 

^IF - V 2 # e = 0, (4) 

where a is the conductivity and /i is the permeability [3]. One should complement Eq. (JJJ) 
with the Faraday condition around the wall 

d x H e {0 + , y, t) - d x H e {0~, y, t) = 2al s v x (t), (5) 

where v x (t) is the domain wall velocity and we have assumed c = oo. Eqs. (J3J and (J3J are 
solved with the boundary condition H e = at the sample surface. The mean pressure on 
the wall is obtained as 

„ 21. r b / 2 



J_ b2 H e (0,y,t)dy = J dt'f(t-t')v x (t'). (6) 



b J -b/2 



A direct calculation of the response function, similar to Ref. [23], leads to 




(7) 



where 9 2 is the Jacobi Elliptic Function and the relaxation times are r n = fiab 2 /[(2n + l)vr] 2 . 

The equation of motion for the domain wall is obtained by equating the eddy current pres- 
sure P e to the pressure P a = 2I S H exerted by the applied field corrected by demagnetizing 
effects and pinning. The left-hand side of Eq. (J2J) is recovered by considering the leading con- 
tributions at short and long times, involving the largest relaxation time r = r = fxab 2 /tt 2 . 
The constants can be estimated as 



In the frequency domain, we can formally rewrite Eq. © as P e = {(3 + iu>M*)v x , with 
the effective mass, in the low frequency limit, given by 





(8) 



with a « 0.886 and 7 = £~ l/(2n + l) 3 « 1.05. 



M* 




2 7 



(9) 
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FIG. 1: Schematic picture of the experimental setup. A ribbon of dimensions a x b x 
c is placed in a solenoid (not shown) yielding a magnetic field in the vertical direction. The 
magnetization reversal process occurs via the displacement of domain walls, such as the one depicted 
here. The arrows indicate the directions of the magnetization. A pickup coil, wound around the 
sample, records a voltage signal v that is composed by distinct pulses, as the one depicted in the 
right panel. 



10 



1.5 



5 
A 
P 

V 



0.5 







1 

- 


1 1 


1 1 1 1 

- 


/t 
/V 
/'/ 
/7 




VX 


/■"/ 

- f / 
II 




\\ - 


ll 

$1 
III 


T=100ns 






T 2()(),us 


\ 




— - T=300ns 




- 


— T=400ns 




1 


1 , 





0.2 0.4 0.6 0.8 

t/T 



1.5 



0.5 







1 1 

f' 

Jf 
Jf 


1 1 1 1 


1 1 1 


'/ 
I 

3 1 

// 

// / 
f '/ 

hi 
ul 
in 




\\ 
^\ 

v.\\ 
v.\\ 


— 1/ 
i! 


— T/x=5 




ul 
ul 


-— T/x=7.5 




- 


T/T- 10 




el 


— - T/x=12.5 




' , i 


1 , 1 


1 



1 0.2 0.4 0.6 

t/T 



0.8 



FIG. 2: Comparison of the average pulse shapes in the experiments and in the model. 

On the left panel we report the average pulse shape obtained from Barkhausen noise measurements 
in a partially crystallized Fe64Co2iBi5 ribbon. The shapes for different durations T are normalized 
and rescale quite well, apart from a small systematic variation related to the asymmetry. On 
the right panel, we report the corresponding shapes obtained from the model. The normalization 
constant is given by N = Jq dt(v(t,T)}/T. 
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FIG. 3: Skewness of avalanches in experiments. The skewness of the pulse shape as a 
function of the pulse duration obtained from experiments. The data display a peak at Tp Xp ~ 200/Us 
indicated by an arrow. 
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FIG. 4: Skewness of avalanches in the model. The skewness of the pulse shape obtained in 



the domain wall model for different values of eddy current time constant r is plotted as a function 



of the pulse duration T. We also report a similar measurement in the case of conventional positive 



domain wall mass. In the latter case the skewness is negative (rightward asymmetry). 
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